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ABSTRACT 

We present computations of the ionization structure and metal-absorption properties of thermally 
conductive interface layers that surround evaporating warm spherical clouds, embedded in a hot 
medium. We rely on the analytical steady-state formalism of Dalton and Balbus to calculate the 
temperature profile in the evaporating gas, and we explicitly solve the time-dependent ionization 
equations for H, He, C, N, O, Si, and S in the conductive interface. We include photoionization by 
an external field. We estimate how departures from equilibrium ionization affect the resonance-line 
cooling efficiencies in the evaporating gas, and determine the conditions for which radiative losses 
may be neglected in the solution for the evaporation dynamics and temperature profile. Our results 
indicate that non-equilibrium cooling significantly increases the value of the saturation parameter 
Co at which radiative losses begin to affect the flow dynamics. As applications we calculate the ion 
fractions and projected column densities arising in the evaporating layers surrounding dwarf-galaxy- 
scale objects that are also photoionized by metagalactic radiation. We compare our results to the UV 
metal-absorption column densities observed in local highly-ionized metal-absorbers, located in the 
Galactic corona or intergalactic medium. Conductive interfaces significantly enhance the formation 
of high-ions such as C 3+ , N 4+ , and 5+ relative to purely photoionized clouds, especially for clouds 
embedded in a high-pressure corona. However, the enhanced columns are still too low to account for 
the O VI columns (~ 10 14 cm" 2 ) observed in the local high- velocity metal- ion absorbers. We find that 
column densities larger than ~ 10 13 cm" 2 cannot be produced in evaporating clouds. Our results do 
support the conclusion of Savage and Lehner, that absorption due to evaporating O VI likely occurs 
in the local interstellar medium, with characteristic columns of ~ 10 13 cm" 2 . 

Subject headings: ISM:general - atomic processes - conduction - plasmas - quasars:absorption lines 



1. INTRODUCTION 

The physics of conductively evaporating clouds has 
applications to a variety of astrophysical environments. 
This includes the study of the Local Cloud and Lo- 
cal Bubble (Slavin 1989; Smith & Cox 2001, Slavin & 
Frisch 2002; Savage & Lehner 2006; Jenkins 2009), cloud 
evaporation in the interstellar medium (e.g. Bohringer 
& Hartquist 1987; Dalton & Balbus 1993; Nagashima et 
al. 2007; Slavin 2007; Vieser & Hensler 2007), supernovae 
remnants (e.g. Shelton 1998; Slavin & Cox 1992; Balsara 
et al. 2008), active galactic nuclei (e.g. McKee & Begel- 
man 1990), and the galaxy formation process (Nipoti & 
Binney 2007). 

In this paper, we present new computations of the non- 
equilibrium ionization states and metal absorption line 
signatures of thermally conductive interfaces that sur- 
round evaporating warm (~ 10 4 K) spherical gas clouds 
embedded in a hot ionized medium (HIM, > 10 6 K). In 
such interfaces, the warm clouds undergo steady evapora- 
tion, while heat from the hot ambient medium flows into 
the clouds (Cowie & McKee 1977; McKee & Cowie 1977). 
We compute the integrated metal-ion column densities 
through the conductive interfaces, for comparison to ab- 
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sorption line observations. We consider the effect of 
photoionization by the metagalactic radiation field on 
the non-equilibrium ionization states in the conduction 
fronts. 

As matter flows from the warm cloud toward the hot 
ambient medium, its ionization state changes continu- 
ously. Non-equilibrium effects become significant when 
the ionization time-scale is long compared to the rate of 
temperature change. The gas then tends to remain "un- 
derionized" compared to gas in ionization equilibrium at 
the same temperature. Because the energy losses in the 
evaporating gas are dominated by atomic and ionic line 
emissions from many species, the cooling efficiencies are 
affected by the non-equilibrium abundances. For under- 
ionized gas, cooling is enhanced. 

Our theoretical work is motivated by recent ultravio- 
let and X-ray absorption line spectroscopy of hot gas in 
the Galactic halo, around higher redshift galaxies, and 
in intergalactic environments (e.g. Sembach et al. 2003; 
Collins et al. 2005; Savage et al. 2005; Tumlinson et 
al. 2005; Fang et al. 2006; Stocke et al. 2006; Narayanan 
et al. 2010). We are also motivated by observations 
of discrete ionized "high-velocity clouds" (Sembach et 
al. 1999, 2000, 2002, 2003; Murphy et al. 2000; Wakker 
et al. 2003; Collins et al. 2003, 2004, 2007; Foxet al. 2005, 
2006). These observations indicate that the ionized gas 
is not in equilibrium ionization (Sembach et al. 2002; 
Fox et al. 2004, 2005; Gnat & Sternberg 2004), and non- 
equilibrium processes such as time-dependent radiative 
cooling (e.g. Gnat & Sternberg 2007; Sutherland & Do- 
pita 1993), shock ionization (e.g. Allen et al. 2008; Gnat 
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& Sternberg 2009), or ionization in conductive interfaces 
(this paper) may be at play in these objects. 

Theoretical models of cloud evaporation in a hot 
medium have been studied extensively. Cowie & Mc- 
Kee (1977) developed an analytical solution for the mass 
loss rate and temperature profile in static conductive 
interfaces around spherical clouds. In their solution, 
they included the possible effects of heat flow satura- 
tion that occurs when the electron mean free path is 
not small compared to the temperature scale height (see 
§2). The saturation in the flow can be measured by 
one global parameter, the "saturation parameter" <jq. 
McKee & Cowie (1977) studied the effects of radiative 
losses on the interface. In particular, they demarcated 
the range of saturation parameters for which radiative 
losses may be neglected. The situation becomes much 
more complicated if the dynamics of the evaporation are 
affected by magnetic fields (e.g. Slavin 1989; Borkowski, 
Balbus & Fristrom 1990) or if the evaporation is time- 
dependent (e.g. Shelton 1998; Smith & Cox 2001; Vieser 
& Hensler 2007). 

Non-equilibrium ionization in conductive interfaces 
was considered previously in several studies (Ballet, Ar- 
naud & Rothenflug 1986; Bohringcr & Hartquist 1987; 
Dalton & Balbus 1993; Slavin 1989; Slavin & Cox 1992; 
Shelton 1998; Smith & Cox 2001). However, most of 
these works explored only a limited range of parameters 
(with -Phim/^b ~ 10 4 cm~ 3 K and -R c i ou d ~ 1 pc), ap- 
propriate for interstellar clouds. None of the previous 
computations included photoionization by an external 
radiation field. 

Here we calculate the ion fractions created in the con- 
ductive interfaces surrounding evaporating gas clouds. 
We consider pressures between 0.1 and 10 4 cm -3 K, cloud 
radii between 1 pc and 100 kpc, and ambient tempera- 
tures between 5 x 10 5 K and 10 7 K. We study how the 
metal absorption line properties depend on the interface 
parameters. We consider the effects of non-equilibrium 
cooling and reexamine the criterion derived by McKee 
& Cowie (1977) for the neglect of radiative losses in the 
evaporative flows. 

As applications, we then focus on the example of the lo- 
cal "ionized high- velocity clouds" . Over the past decade, 
UV absorption line observations have revealed a popula- 
tion of local high- velocity metal-ion absorbers (Sembach 
et al. 1999, 2000, 2002, 2003; Murphy et al. 2000; Wakker 
et al. 2003; Collins et al. 2004; Fox et al. 2005, 2006). 
The observations, carried out with the Goddard High 
Resolution Spectrometer (GHRS) and Space Telescope 
Imaging Spectrograph (STIS) on board HST, and more 
recently with the Far Ultraviolet Spectroscopic Explorer 
(FUSE), indicate the presence of many high- ionization 
metal absorption lines, including Si III A 1206.5, Si IV 
AA 1393.8, 1402.8, and C IV AA 1548.2, 1550.8 with lit- 
tle or no associated H I (N H i < 10 15 - 10 16 cm" 2 ). Of 
particular interest is O VI A 1031.9, for which absorbing 
column densities of order 10 14 cm -2 have been measured 
(Collins et al. 2004; Fox et al. 2006; Sembach et al. 2000, 
2003). 

Many of these ionized absorbers can be identified with 
the diffuse nearby H I high-velocity cloud complexes, 
such as the Magellanic Stream or Complex C (Sem- 
bach et al. 2003). However, some of these absorbers do 
not appear directly linked to such structures, and could 



be clouds at much larger distances, perhaps pervading 
the Local Group. The absorbers towards Mrk 509 and 
PKS 2155-304 are examples of such isolated high- velocity 
metal-ion absorbers (Sembach et al. 1999; Collins et 
al. 2004; Fox et al. 2005). As a guide to these observa- 
tions, we list in Table [1] the data presented by Collins et 
al. (2004) for the metal-ion absorbers towards Mrk 509. 

In Gnat & Sternberg (2004, hereafter GS04), we mod- 
eled the high-velocity metal ion absorbers as photoion- 
ized gas clouds associated with low-mass dark matter 
"minihalos" (see also Kepner et al. 1999). For this pur- 
pose, we studied the metal photoionization properties of 
hydrostatic gas clouds embedded in gravitationally dom- 
inant dark-matter halos, that are photoionized by the 
present-day cosmological metagalactic radiation field, 
and are pressure-confined by an external hot medium. 
We considered minihalo models for dwarf-galaxy-scale 
objects and for the lower mass compact H I high veloc- 
ity clouds (CHVCs), based on the properties derived by 
Sternberg, McKee, k WolfircQ (2002; see also Giovanelli 
et al. 2010). 

For low mass objects (M V - 1T ~ 10 8 M Q ) embedded in 
the relatively high-pressure environment of the Galactic 
corona (P ~ 50 cm~ 3 K), we found that the ionization 
parameter is too low to efficiently produce C IV and other 
high ions. However, for the more massive dwarf-galaxy- 
scale halos (M vir ~ 2 x 10 9 M Q , P < 1 cm -3 K) em- 
bedded in the low-pressure intergalactic medium (IGM) , 
we found that their ionized envelopes are natural sites 
for the formation of high-ions. Photoionized envelopes 
of dwarf-galaxy-scale halos could be detectable as UV 
metal-line absorbers, with ionization states similar to 
those observed in the ionized high-velocity absorbers. 
O VI remains an important exception, as it is inefficiently 
produced by photoionization. 

The strong O VI absorption in the high- velocity metal- 
ion absorbers strongly suggests that photoionization is 
not the only ionization mechanism at play in these ob- 
jects (Sembach et al. 1999; Fox et al. 2004, 2005; GS04). 
The contribution of additional ionization processes may 
produce high-ions that are suppressed in the purely pho- 
toionized clouds. 

The fact that minihalo clouds may be embedded in a 
hot plasma, such as the Galactic corona or the intergalac- 
tic medium at various redshifts, suggests the possibil- 
ity of significant collisional interactions between the hot 
gas and the warm clouds. Here we use our evaporating 
cloud models to study the interface layers surrounding 
dwarf-galaxy-scale minihalos (0.1 < P/k-g, < 50 cm~ 3 K, 
-Rcioud ~ kpc, Thim > 10 6 K, Z = 0.1 solar) for compar- 
ison with the observations of the high- velocity metal ion 
absorbers. As we describe below, we include the effect 
of photoionization by the metagalactic external radiation 
on the non-equilibrium ionization properties of the gas in 
the interface. 

In §2 we describe the basic equations and numeri- 
cal method. In §3 we present a set of computations 

5 Sternberg, McKee, & Wolfire (2002) constructed explicit mini- 
halo models for Local Group dwarf galaxies, based on the observed 
H I properties of Leo A and Sag DIG. They then searched for 
minihalo models for the H I CHVCs, that would simultaneously 
account for the range of observed H I column densities, the exis- 
tence of multiphased (cold/warm) cores, and the total number of 
objects. 
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TABLE 1 

Column-Density Measurements in MrK 509 



Ion 


-300 km s- 1 


-240 km s- 1 


CII 
N I 
Si II 
S II 


4.79t°:gS x 10" 

< 3.39 x 10 13 

< 1.23 x 10 12 

< 1.00 x 10 14 


< 2.09 X 10 13 

< 2.40 X 10 13 
8.5lt 2 ;°5 x 10 12 

< 7.08 x 10 13 


civ 

N V 
O VI 
Si III 
Si IV 
S III 


1A1±°H x 10 14 

< 1.74 x 10 13 
8.5ltg;g x 10" 

2 - 04+ oil x 1Ql3 
3.39±i;il x 10 13 

< 8.51 x 10 13 


3.39±°;« x 10 13 

< 1.20 x 10 13 
7.76+g-g x 10 13 
2.75iSiS x 10 12 

< 3.09 x 10 12 

< 6.03 x 10 13 



Note. — Column densities for two velocity components toward 
Mrk 509 by Collins et al. (2004). The error estimates are the la- 
estimates, and the upper limits are 3<r limits. 

of the ion fractions in conductive interfaces, and we 
discuss the resulting metal absorption column densities 
for the dwarf-galaxy-scale models in §4. We find that 
conductive interfaces significantly enhance the forma- 
tion of high ions such as C 3+ , N 4+ and 5+ relative to 
purely photoionized clouds, especially for clouds embed- 
ded in a high-pressure medium. However, the enhanced 
columns are still too low to account for the O VI columns 
(~ 10 14 cm -2 ) observed in the high-velocity metal-ion 
absorbers. We find that O VI column densities larger 
than ~ 10 13 cm~ 2 cannot be produced in evaporating 
clouds. Our models do support the conclusion by Sav- 
age & Lehner (2006) that evaporating O VI absorption 
occurs in the local ISM, with characteristic columns of 
~ 10 13 cm -2 . We summarize in §5. 

2. BASIC EQUATIONS AND PROCESSES 

2.1. Evaporation 

We are interested in studying the evolving non- 
equilibrium ionization states in thermally conductive gas 
that evaporates from warm purely photoionized clouds 
into a hot ambient medium. In our treatment of the ion- 
ization in the evaporating layers we use an analytic solu- 
tion for the temperature profiles in the conductive inter- 
faces surrounding steady-state, spherical, non-magnetic 
clouds. The temperature profiles are obtained by solving 
the equation of energy conservation in the flow, balanc- 
ing the outward energy flux with the inward heat flux. 
Within the interface, the heat flux changes from a clas- 
sical diffusive form, 



q cl = -kVT 



(1) 



(Spitzer 1962) that applies when the electron mean free 
path is small compared to the temperature scale height 
(T/VT), to a saturated form, 



|q S at| = h(j) s pc 



(2) 



(Cowie & McKee 1977) when the mean free path is com- 
parable to, or greater than, the temperature scale height. 
In the above expressions, n is the thermal conductivity^, 

6 k = 1.84 X 10- 5 T 5 / 2 (lnA)- 1 erg s" 1 cm" 1 K" 1 where In A = 
29.7 + ln(nr 1/2 T 6 ), n e is the electron density, and T 6 = T/(10 6 K) 



p is the mass density, c is the sound speed, and <p s is a 
factor of order unity that allows for various uncertainties 
in the estimation of q sa ± . 

The local saturation in the flow is defined by the ratio 
(Cowie & McKee 1977), 

o = , 

<7sat 



(3) 



and is proportional to the ratio of the mean free path to 
the temperature scale height. 
For a classical diffusive interface, the rate of mass loss 



m c i 



16717ii?K H IM 

25A^ 



(4) 



(Cowie & McKee 1977), where R is the cloud radius, 
p is the mean mass per particle, fcs is the Boltzmann 
constant, and subscript "HIM" refers to the surrounding 
hot ionized medium. As the flow saturates, the heat flux 
is limited to a maximal rate given by equation (2), and 
the classical mass loss rate is reduced by a factor uj such 
that m = ujrh c \. The global level of saturation in the flow 
is parametrized by the saturation parameter, 



cto 



2khimThim 



oc 



T 2 



1 HIM 



25<f) s 



(5) 



(Cowie & McKee 1977). In this expressions T is the 
temperature, and P — pk^T/ p is the gas pressure. 

To obtain the steady state temperature profile in the 
conductive interface, the equation of energy conserva- 
tion, 



f(v 2 + 5c 2 ) 



+ (ti 2 A - nT) = - V • q (6) 



must be solved, where A is the cooling efficiency (erg cm 3 
s _1 ), and T is the heating rate (erg s _1 ). The first term 
on the left hand side of equation © is a bulk kinetic 
energy term; the second is due to internal energy and 
PdV work; and the next two terms account for radiative 
losses (cooling) and gains (heating) . This outward energy 
flux is balanced by an inward heat flux represented by the 
right hand side of equation ((6|. 

We rely on the formalism developed by Dalton & Bal- 
bus (1993, hereafter DB93) for determining the tempera- 
ture profile in the interface. The DB93 formalism allows 
the thermal conductivity to change continuously from 
the classical diffusive form to a saturated form, by using 
an effective heat flux q, defined as the harmonic mean of 
the local classical and saturated fluxes, 



(7) 



1 _ J_ 1 

Q ~ Vol 9sat 

DB93 then solve for the steady state temperature profile 
that develops in the interface layer between the warm 
cloud and the hot medium surrounding it. 

When solving the energy equation, DB93 neglect the 
kinetic energy term, as well as the radiative losses (cool- 
ing) term. We discuss both assumptions in detail below. 
With these simplifying assumptions, the equation of en- 
ergy conservation becomes, 



V • ( -pvc* 



-Vq. 



(8) 



(Spitzer 1962). 
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The temperature profile in the conductive interface^ then 
depends only on the global saturation parameter, do, 
which is a function of the cloud radius, and the ambi- 
ent pressure and temperature, as given by equation ([5])- 
Temperature profiles for various values of Uq are shown 
in Figured] 




Fig. 1. — Scaled temperature, T/Tjiim versus scaled distance 
r/R, computed using the DB93 analytical solutions, for a range of 
saturation parameters. 



2.1.1. On the Validity of the Subsomc Approximation 

The kinetic energy term may be neglected when v 2 <C 
5c 2 (see equation [5]), or equivalently when M 2 <C 5. If 
the kinetic energy is small relative to the internal energy, 
the flow is isobaric, with a constant gas pressure P = 

Phim^bThim/ 

In classical flows, M is always less than unity, and 
the subsonic approximation holds. However, as the flow 
saturates, the maximal Mach number increases. In fully 
saturated flows, the heat flux is given by equation (2), 
and equation (6) reads 



(9) 



(10) 



5/2 /wrV(l + M 2 /5) - 50 s/ oeV = 0, 

This implies that in fully saturated flows 

M(l + M 2 /5) = 2<j) s 

(e.g. Cowie & McKee 1977). Therefore, for <p s < 0.6, the 
subsonic approximation holds throughout the flow. 

As a consistency check, we also consider the local Mach 
number in the analytic profiles derived by Dalton & Bal- 
bus neglecting the kinetic energy term. In this case 



M _ T i/2 k B u f l§TTy} /2 K. mM R 
P Anr 2 V 



(id 



(DB93 equation B3). Given a set of temperature pro- 
files appropriate for a range of saturation parameters (see 

7 DB93 write the solution for the scaled radius, y = r/R, as a 
function of the scaled temperature, r = T/Tjjimi as j/ = j^uix'r 1 / 4 ' 

^5/ll(x)I-(6/ll)(XT 11/4 ) + g 6 /ll(xT 11/4 )-r-( 6 /ll)(x) j ,. 

K 5 /ll(X- 11 / 4 )/5/ll(x)-^5/ll(x)/5/ll(xr 11 / 4 ) • lntmSeX P reS - 

sion x = (jx)(^) 1//2 ' ^ are * ne modified Bessel functions of the 
first kind, and K are the modified Bessel functions of the second 
kind. 



Figure Q} , this equation shows that the maximal Mach 
number is an increasing function of <tq. With the ki- 
netic energy term neglected, the maximal Mach number 
approaches 2(f> s in highly saturated flowt0, implying that 
for 4> s < 0.5 the subsonic approximation holds even for 
large a . 

Experimental evidence suggests that <f> s ~ 0.3 (see Bal- 
bus & McKee 1982). Throughout this paper we therefore 
assume <f> s — 0.3, consistent with the use of subsonic ap- 
proximation. 

2.1.2. On the Validity of the Non-Radiative Approximation 

The fate of warm clouds that are embedded in a hot 
ambient medium depends of the delicate balance between 
the conductive heating and the radiative losses in the 
interface. Evaporation takes place if the radiative losses 
are small compared to conductive heating. If, on the 
other hand, cooling overcomes the conductive heating, 
the evaporation stops, and instead the ambient medium 
condenses onto the warm cloud. 

McKee & Cowie (1977) used an approximate analytical 
method to determine the critical cloud radii at which ra- 
diative losses balance the conductive heating. They used 
fits to the equilibrium cooling functions of solar metal- 
licity gas to derive these critical radii as a function of 
the external temperature and density. McKee & Cowie 
found that radiative losses are important only at low sat- 
urations, (To < 0.03/0,,. Such small values of <jq are gen- 
erally out of the regime "occupied" by dwarf-galaxy scale 
halos (for which cr > 0.2). 

Our models differ from those presented by McKee & 
Cowie. First, we explicitly follow the non-equilibrium 
ion fractions in the interface. As mentioned above, the 
outflowing gas may be underionized relative to gas in 
equilibrium at the same temperature. Underionized gas 
tends to radiate more efficiently (McCray 1987; Gnat & 
Sternberg 2007), thus increasing the contributions of ra- 
diative losses. Second, we include photoionization by an 
external background radiation field. This both increases 
the ionization level in the gas and provides heating, both 
effects acting to reduce the importance of radiative losses. 
Finally, we assume gas metallicity of 0.1 solar, thus re- 
ducing the cooling efficiency. The relative importance of 
these three effects must be estimated numerically. 

Using our numerical results we check, for each model, 
whether the cmissivity associated with the heat flux is in- 
deed everywhere larger than the radiative cooling. Given 
our non-equilibrium ion fractions (including photoioniza- 
tion, see section |2~2"|) . we use the cooling functions in- 
cluded in Cloudy (version 07.02, Ferland et al. 1998) to 
compute the local cooling efficiencies. We then compare 
these values with V • q. Figure [2] displays our results. 

8 DB93 argued, incorrectly, that the maximal Mach number is a 
decreasing function of <jq . They used a specific set of physical vari- 
ables (P/k B = 10 4 cm" 3 K, T H im = 7 x 10 5 K, i?, = 1 pc) to derive 
the numerical expression M = 0.66a; \JT /Thim / (r / R) 2 (equation 
B5 in DB93). They then used this expression to evaluate the Mach 
number for a range of saturation parameters, and concluded that 
the maximal Mach number is a decreasing function of (To- How- 
ever, the numerical coefficient in their expression was computed for 
interface variables implying ctq =0.32 and is therefore not valid for 
other saturation parameters. 

9 As opposed to max(M) = 1.42 (1.0, 0.56) resulting from 
equation (10) which includes the kinetic energy term, for <j>s = 1.0 
(0.6, 0.3). 
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Fig. 2. — (a) Maximum of n 2 A/V ■ q versus saturation param- 
eter cr . The black (P/k B < 10 3 cm -3 K) and red (P/k B = 
10 4 cm -3 K) points are for Z = 0.1 solar metallicity gas. For 
P/^B = 10 4 cm -3 K the neutral fraction on the cloud surface is 
considerable (xjji = 0.574). The gray circles are for solar metallic- 
ity gas (all pressures), (b) Fraction of interface temperature range 
over which the cooling is locally larger than V ■ q. Triangles indicate 
models for which this fraction is smaller than 0.003. 

Panel (a) shows the maximal ratio of the radiative to 
conductive emissivity in each mode£3 as a function of 
the saturation. Panel (b) shows the fraction of the in- 
terface temperature over which the cooling is larger than 
the conductive heating as a function of saturation. The 
points (black and red) are for our assumed heavy element 
abundances of Z = 0.1 times solar. The gray circles show 
results for solar metallicity gas, for a direct comparison 
with the critical <jq derived by McKee & Cowie (1977). 

For Z = 0.1 times solar, both panels show that for 
(To < 0.2 (for 4> s = 0.3) radiative losses are significant 
and may not be neglected. For larger values of <r the 
importance of radiative cooling is generally smaller, and 
depends on the specific parameters of the evaporating 
cloud. 

For <7o ^ 0.2, enhanced cooling and larger values of 
n 2 A/V • q are associated with interfaces of higher gas 
densities, with larger neutral fractions at the cloud sur- 
face. We find that for our assumed 10 4 K-photoionized 
clouds, the neutral fraction on the surface is significant 

10 We only display the maximum for T < 0.9Thim- As the 
temperature approaches the ambient temperature, the heat flux 
approaches zero but the cooling does not. 



when P > 10 cm~ 3 K. For example, for pressures of 
P = 10 4 , 10 3 , and 100 cm~ 3 K, the neutral fractions are 
0.57, 0.20, and 0.03, respectively. Neutral gas cools very 
efficiently, and in underionized flows the lingering and 
enhanced contribution of H I Lya cooling significantly 
affects the energy balance. 

To study the impact of the surface ionization on the 
flow dynamics, we distinguish between "neutral cloud" 
models for which P/fes ^ 10 4 cm~ 3 K (red points), 
and "ionized cloud" models with P/k-Q < 10 3 cm~ 3 K 
(black pointsjE]. Figure [2] shows that for <jq > 0.2, 
radiative losses are negligible for bounding pressures 
< 10 3 cm -3 K. However, in more neutral clouds cool- 
ing may be significant even at high saturations. 

The fact that the limiting saturation parameter that 
we find is similar to the one derived by McKee & Cowie 
(1977, do = 0. 03/(f> s ) - despite the 10 times lower metal- 
licity that we use, and the inclusion of photoioniza- 
tion - implies that the cooling efficiencies in our non- 
equilibrium interfaces are considerably larger than their 
equilibrium counterparts. In fact, we find that Lya from 
underionized He + dramatically increases the cooling ef- 
ficiency in the interface. We discuss this point in detail 
in section [3l 

For Z = 1 solar, we find that radiative losses are im- 
portant up to a critical saturation parameter of do ~ 0.5 
(for <f> s = 0.3), and become less significant for larger val- 
ues of (To • Since underionized Hc + is one of the dominant 
coolants in the interface, increasing the metallicity from 
a tenth-solar to solar changes the critical saturation pa- 
rameter by a factor smaller than 10. Figure [5] shows that 
this factor is ~ 3. 

For the tenth solar metallicity dwarf-galaxy-scale mod- 
els that we use for comparison with the observations of 
the high- velocity metal-ion absorbers, the radiative losses 
are, in fact, everywhere smaller than the conductive heat- 
ing. Some examples are shown in Figure [31 in which we 
compare V • q (solid curves) with the radiative losses 
(dashed curves) in three interface models. The left hand 
panels are for interfaces embedded in a high pressure 
corona, with P/k-Q = 50 cm -3 K. We show results for a 
small, marginally saturated cloud (R = 1.3 kpc, Co ~ 3) 
in panel (a), and for a classical interface surrounding a 
larger cloud (R = 7.5 kpc, ern ~ 0.6) in panel (b). In 
panel (c) we present a highly saturated interface, with 
do ~ 270, surrounding a 7.5 kpc cloud, embedded in a 
low pressure intergalactic medium (P/kB = 0.1 cm -3 K). 
In all cases, the radiative losses are everywhere smaller 
than the emissivity associated with the heat flux. These 
plots confirm that the significance of radiative losses 
is larger for lower saturation parameters, and becomes 
small at larger ctq. 

2.2. Ionization 

Given the DB93 temperature profile, we solve for the 
non-equilibrium abundances of the various species in the 
interface. We follow the time-dependent ionization in the 
evaporating gas, given the initial photoionization equi- 
librium state in the warm cloud. We consider all ioniza- 
tion stages of the elements H, He, C, N, O, Si, and S. 

11 P/k-g, = 10 3 cm -3 K corresponds to an ionization parameter 
of 10 on the cloud surface (see Section 3, equation [TH}. 
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Fig. 3. — Comparison of the radiative losses and V ■ q (erg s 1 cm 3 ), for three interface models, with <ro of 3.34 (a), 0.58 (b), and 267 
(c). 



The temperature-dependent ionization and recombina- 
tion processes that we include are collisional ionization 
by thermal electrons (Voronov 1997), radiative recombi- 
nation (Aldrovandi & Pequignot 1973; Shull & van Steen- 
berg 1982; Landini & Monsignori Fossi 1990; Landini & 
Fossi 1991; Pequignot, Petitjean, & Boisson 1991; Ar- 
naud & Raymond 1992; Verner et al. 1996), dielectronic 
recombination (Aldrovandi & Pequignot 1973; Arnaud & 
Raymond 1992; Badnell et al. 2003, Badnell 2006; Col- 
gan et al. 2003, 2004, 2005; Zatsarinny et al. 2003, 2004a, 
2004b, 2005a, 2005b, 2006; Altun et al. 2004, 2005, 2006; 
Mitnik & Badnell 2004), and neutralization and ioniza- 
tion by charge transfer reactions with hydrogen and he- 
lium atoms and ions (Kingdon & Ferland fitfi. based on 
Kingdon & Ferland 1996, Ferland et al. 1997, Clarke et 
al. 1998, Standi et al. 1998; Arnaud & Rothenflug 1985). 

The time-dependent equations for the ion abundance 
fractions, Xi, of element m in ionization stage i are, 



v dxi/dr = Xi 



[Qi-irie + rj_i + fc] i i _ 1 n H + 
+^ < L 1 n Ho +] 
+x i+ i [a l+1 n c + kf i+1 n s a + kf? +1 n He o 
-X{ [(qi + ai)n c + + fe^n B 



(12) 



In this expression, v is the gas velocity at radius r in the 
evaporating flow. The parameters qi and a% are the rate 
coefficients for collisional ionization and recombination 
(radiative plus dielectronic), and fc^, fcj*, fcjj 3 , and fcj?^ 
are the rate coefficients for charge transfer reactions with 
hydrogen and helium that lead to ionization or neutral- 
ization. The quantities n^o , , hh c o , «h + > and n e 
are the particle densities (cm~ 3 ) for neutral hydrogen, 
ionized hydrogen, neutral helium, singly ionized helium, 
and electrons, respectively. Ti are the photoionization 
rates of ions i, due to externally incident radiation. 

For the external radiation, we use the Sternberg, Mc- 
Kee & Wolfire (2002) fit for the present-day metagalactic 

12 See: |http://www-cfadc.phy.ornl.gov/astro/ps/data/cx| /hy- 
drogen/rates/ct.html 
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v ° 
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(13) 

to compute Ti. In this expression J„o — 2 x 10 23 erg 
s^ 1 cm -2 Hz" 1 sr" 1 , and v$ is the Lyman limit fre- 
quency. This fit, which we plot in Figure HI is based on 
observational constraints for the optical/UV and X-rays, 
(Bernstein et al. 2002; Martin et al. 1991; Chen et al. 
1997) and on theoretical models (Haardt & Madau 1996) 
for the unobservable radiation from the Lyman limit to 
- 0.25 keV. 

Recent observations and modeling have suggested that 
the spectral slope of the metagalactic radiation may be 
shallower near the Lyman limit, with J v varying as j/~ 1,8 , 
and steepen gradually with increasing energy (e.g. Shull 
et al. 1999; Telfer et al. 2002; Scott et al. 2004; Shull et 
al. 2004; Zheng et al. 2004; Faucher-Giguere et al. 2009). 
In Appendix A we consider a radiation field (also plotted 



in Figure H]) that varies as 



-1.725 



from the Lyman limit 



to 0.25 keV, and investigate the sensitivity of our results 
to the spectral slope. For reference, we also plot in Figure 
4 the Faucher-Giguere et al. (2009) model computation 
of the present-day metagalactic field. In Appendix A we 
show that the column densities of the high ions (O VI and 
N V) are insensitive to the spectral slope at the Lyman 
limit, while the column densities of lower ions (e.g. Si III) 
may be affected at a level of < 0.5 dex, mostly for low 
external bounding pressures. 

The present-day metagalactic radiation field. The solid 
dark line is for the standard field as represented by Stern- 
berg et al. (2002; see text and Equation fT5|) . the gray 
line is for the Faucher-Giguere et al. (2009) field, and 
the dashed line is for the shallower "maximal field" con- 
sidered in Appendix A (Gnat & Sternberg 2004). 

For each element to, the ion fractions Xi = 
ni im /(nnA m ) satisfy 



E 



Xi 



1 



(14) 
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Fig. 4. — The present-day metagalactic radiation field. The solid 
dark curve is for our standard fiel d as represented by Sternberg et 
al. (2002; see text and Equation |13> . the gray curve is for the 
Faucher-Giguere et al. (2009) field, and the dashed curve is for 
the shallower "maximal field" considered in Appendix A (Gnat & 
Sternberg 2004). 



where ni jjn is the density (cm -3 ) of ions in ionization 
stage i of element m, «h is total hydrogen density, and 
A m is the abundance of element m relative to hydrogen. 
The sum is over all ionization stages of the element. 

For the atomic elements that we include, equa- 
tions (12) are a set of 54 coupled ordinary differential 
equations (ODEs). For the initial conditions, we assume 
that the cloud surface is in photoionization equilibrium 
with the metagalactic filed, and has a temperature of 
10 4 K. We advance the solutions in small spatial steps 
Ar, where r is the current distance from the cloud cen- 
ter (r > R), which is associated with a temperature T(r) 
through the analytical solutions outlined by DB93. We 
integrate equations (12) over the interval Ar using the 
same Livermore ODE solveiEl (Hindmarsh 1983) we used 
in Gnat & Sternberg (2007). We verified that our code 
converges to the equilibrium abundances given a constant 
temperature, by comparing our solutions to those found 
by solving the set of algebraic equations for X{ obtained 
by setting d,Xi/dr = in equation (12), as appropriate 
for steady state. We verified our time-dependent com- 
putations by comparing with previous non-equilibrium 
models (e.g. Sutherland & Dopita 1993), as discussed in 
detail in Gnat & Sternberg (2007). 

3. METAL IONS FRACTIONS AND COOLING IN 
EVAPORATING CLOUDS 

We computed the ionization states of the elements H, 
He, C, N, O, Si, and S in thermally conductive inter- 
faces surrounding evaporating warm gas clouds. We cal- 
culate the non-equilibrium ionization states as a func- 
tion of distance from the warm cloud surface. We con- 
sider a grid of models with pressures between 0.1 and 
10 4 cm" 3 K, cloud radii between 1 pc and 100 kpc, and 
ambient (HIM) temperatures between 5 x 10 5 and 10 7 K. 
We assume heavy element abundance of 0.1 times solar, 
as listed in Table [5] We begin with a brief discussion 
of our non-equilibrium conductive interfaces, and later 
focus on models for dwarf-galaxy-scale objects. 



See: http://www.netlib.org/odepack/ 



TABLE 2 
0.1 SOLAR METALLICITY ABUNDANCES 



Element 


Abundance 




(X/H) 


Carbon 


2.45 x 10" 5 


Nitrogen 


6.03 x 10" 6 


Oxygen 


4.57 x 10" 5 


Silicon 


3.24 x 10" 6 


Sulfur 


1.38 x 10"" 6 



As the evaporating gas flows into the external, hot- 
ter, parts of the interface layer, its overall ionization 
state gradually increases. If the gas is heated faster than 
it is ionized, nonequilibrium effects become significant 
and the gas remains "undcrionized" throughout the flow. 
This happens when the heating time-scale, which is set 
by the flow time through the interface, 



ifio 



dr R 3 pmM 

OC ; . 

v m 



(15) 



becomes short compared to the ionization time scale for 
ion i, ti onj i. 

In the absence of external photoionization, the ioniza- 
tion time scale is given by 



t 1 , 



3 — rad 



i 



»eCi(T)' 



(16) 



where Ci(T) is the collisional ionization coefficient for 
ion i at temperature T. The ionization state re- 
mains close to collisional ionization equilibrium when 
iflow/iion,i 3> 1. Nonequilibrium effects become signifi- 
cant when tflowAion.i < 1- With no photoionization, the 
ratio of the flow to ionization time scale, 



-rad 



T 3 / 2 d(T) 



(17) 



is a function just of the saturation parameter and ambi- 
ent temperature. The ion fractions and integrated metal- 
ion column densities then depend on just the two param- 
eters <7o and ?him- 

Here we also include photoionization by the external 
metagalactic radiation field, which increases the level of 
ionization in the gas, and provides heating. The ioniza- 
tion parameter at the surface of a photoionized, 10 K 
cloud is given by, 



4tt 



U = / — dv = 4.3 x 10~ 7 n 

"hc J m hv 



l 

H ' 



(18) 



where v is the Lyman limit frequency, c is the speed 
of light, and nn is the hydrogen density (cm -3 ) at the 
cloud surface. Photoionization is more efficient in lower 
density gas. With photoionization included, the ioniza- 
tion states and integrated column densities depend on 
all three interface parameters, Trim, -P/&B, and R. In 
general, we find that the ionization lag is so great that 
despite the inclusion of photoionization, the gas remains 
highly underionized, not only compared to photoioniza- 
tion equilibrium, but also compared to collisional ioniza- 
tion equilibrium (CIE). 
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Because the energy losses are dominated by atomic and 
ionic emission lines, the nonequilibrium cooling rates also 
differ (are enhanced for underionized gas) compared to 
equilibrium cooling. The lower ionization species that 
remain present at each temperature, most significantly 
He + , are more efficiently excited by the thermal electrons 
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Fig. 5. — Underionized gas in a model with R = 1.3 kpc, 
PaiM/k = 50 cm -3 K, T mM = 2 X 10 6 K. (a) Non-equilibrium 

C 3+ 

ion fraction (dashed curve). The carbon CIE ion fractions are 
shown by the solid curves for comparison, (b) Enhanced cooling 
in the underionized interface. The dashed curve shows the non- 
equilibrium cooling efficiency, and the solid curve shows the CIE 
cooling efficiency (for Z = 0.1 solar). 

We demonstrate the effect of "under-ionization" in 
Figure El for a model with R = 1.3 kpc, Pmu/k — 
50 cm" 3 K, and Trim = 2 x 10 6 K (a = 3.34). In 
the upper panel we show the nonequilibrium C 3+ ion 
fraction versus temperature (dashed curve) , and the CIE 
ion fractions of C 3+ -C 5+ (solid curves, Gnat & Stern- 
berg 2007). The nonequilibrium C 3+ abundance peaks 
at a temperature of ~ 4 x 10 5 K, whereas the CIE C 3+ 
abundance peaks at ~ 10 5 K . For CIE, the dominant 
ion at 4 x 10 5 K is C 4+ . 

The modified ionization states enhance the cooling ef- 
ficiency in this interface. We computed the local cool- 
ing efficiencies using Cloudy (version 07.02, Ferland et 
al. 1998) given our nonequilibrium ion fractions. This is 
shown in panel (b). The CIE cooling efficiency is shown 
by the solid curve (Gnat & Sternberg 2007). The non- 
equilibrium cooling efficiency (dashed curve) is an order 



of magnitude larger in this model, over a significant tem- 
perature range. This is mostly due to He + which persists 
to much higher temperatures in the nonequilibrium gas 
compared to CIE. 
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Fig. 6. — Departures from equilibrium in conductive flows, (a) 
The temperature to which He+ persists, with a fractional abun- 
dance larger than 10%, in the non-equilibrium gas. The CIE value 
is shown by the dashed line. These results are independent of Z. 
(b) Maximal ratio of non-equilibrium to equilibrium cooling versus 
saturation for Z = 0.1 solar (black points), and for Z = 1 so- 
lar (gray circles). Models inside the shaded areas arenot self- 
consistent, as cooling may not be neglected (see section ^. 1.21 1. 



In Figure [5] we demonstrate the effect of non- 
equilibrium ionization for our complete set of models. 
In panel (a) we show the ionization lag for He + . For 
each model, we plot the temperature to which He + per- 
sists, with a fractional abundance larger than 10%, ver- 
sus its saturation parameter. These values are indepen- 
dent of the gas metallicity. The corresponding CIE value 
(~ 10 5 K) is shown by the dashed line. The shaded 
region indicates the range of saturation parameters for 
which no self-consistent evaporating solution exists (for 
Z = 0.1 solar). In this region, radiative cooling may 
not be neglected (see section |2.1.2[) . and results in con- 
densation. For the evaporating models, departures from 
equilibrium ionization can clearly be seen. He + persists 
to higher temperatures than in CIE, despite the contri- 
bution of external photoionization by the metagalactic 
radiation field. For highly saturated models, in which the 
temperature profile is extremely steep, He + may persist 
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to T ~ 10 6 K. 

The impact of these extended ion distributions on the 
cooling efficiencies is shown in panel (b) . For each model 
we plot the maximal ratio of the non-equilibrium- to 
CIE-cooling efficiencies through the interface. The black 
points are for Z = 0.1 times solar, and the gray circles 
are for Z = 1. For the evaporating models, cooling is in- 
deed enhanced by factors ~ 5 to 10 3 relative to the CIE 
cooling efficiency. 

These enhanced cooling efficiencies explain the results 
presented in section 12.1.21 (Figure dj) . We found that for 
solar and 0.1 solar metallicity gas, radiative cooling has 
a small impact on the dynamics of the flow for satura- 
tion parameters cto > 0.5 and cto > 0-2, respectively. 
Our results indicate that non-equilibrium effects signifi- 
cantly increase the value of cto at which radiative effects 
set in. For comparison, in CIE radiative losses in so- 
lar metallicity gas become significant for cto < 0.1 (for 
4> s = 0.3, McKee and Cowie 1977), a factor of ~ 5 lower 
than for nonequilibrium cooling. The reason for this is 
that the cooling efficiencies in the underionized gas are 
larger, due to the enhanced contributions of He + and 
metal-ion resonance-line cooling. 

Does non-equilibrium ionization modify the observa- 
tional signatures of evaporating clouds? In Figure [7f a) 
we compare the equilibrium (dashed curves) and non- 
equilibrium (solid curves) ion distributions versus posi- 
tion in the evaporating layers, for a model with R = 
1.3 kpc, Pmu/k = 50 cm™ 3 K, and Trim = 2 x 10 6 K. 
We show the abundances of C 3+ , N 4 +, 5 +, and Si 3 + 
including the effect of photoionization. The difference 
between the equilibrium and non-equilibrium curves is 
apparent: the non-equilibrium abundances reach their 
peaks further away from the cloud center (at higher tem- 
peratures) as the gas remains underionized, and they 
span considerably larger path-lengths than the equilib- 
rium abundances. These larger path lengths result in 
enhanced column densities. 

Finally, in panel (b) of Figure [7] we demonstrate the 
impact of photoionization by the metagalactic radiation 
field on the ionization states. Here we show a low- 
pressure model with Phim/^ = 0.1 cm™ 3 K, Trim = 
2 x 10 6 K, and R — 7.5 kpc. We compare the ioniza- 
tion state including photoionization by the metagalactic 
radiation field to the ionization in the absence of any 
external radiation. The presence of the external field in- 
creases the level of ionization both in the interface and 
in the surrounding HIM. 

3.1. Metal Ion Fractions in Dwarf-Galaxy-Scale 
Objects 

We now focus on the ion fractions in conductive inter- 
faces surrounding dwarf galaxy-scale halos for compari- 
son with the observed metal-ion absorbers. As a guide 
to the physical properties of dwarf-scale halos, we rely 
on the results presented in GS04, where we considered 
median halos with a virial mass ~ 2 x 10 9 M Q . For the 
dwarf-scale halos, we assume ambient temperatures of 
10 6 K, 2 x 10 6 K, and 10 7 K, gas cloud radii of 3 kpc, 
4.2, and 7.5 kpc, and bounding pressures in the range 
0.1 — 50 cm™ 3 K. For bounding pressures in this range, 
the hydrogen in the 10 4 K photoionized clouds is fully 
ionized at the surface, and the ionization parameter may 



be written as U = 9.7 x 10™ 3 (P/fce)™ 1 . 

In Table [3] we list the input parameters for the dwarf- 
scale models that we consider. These include the ambi- 
ent temperature Thim> the cloud radius R, and the gas 
pressure P/k&. We also list in Table [3] the ionization 
parameter U, saturation parameter cto, the classical (dif- 
fusive) mass loss rate m c i, the mass-loss suppression fac- 
tor uj = m/rhci, the maximal Mach number in the flow, 
max(M), and the evaporation time scale, t ev&p = 1717/171, 
for clouds with warm gas masses of 10 7 M©. All the 
values listed in Table |3] are computed assuming <j) s = 0.3. 

In Figure set 8, we present the ion fractions as a func- 
tion of distance from the cloud surface for the grid of 
dwarf galaxy-scale models. The ambient temperature, 
Thim, cloud radius, R, and bounding pressure, P/fce are 
indicated in each panel. Figure 8 displays the ion frac- 
tions of the low ions C II, N I, O I, Si II, and S II, and the 
high-ions C IV, N V, O VI, Si III, Si IV, and S III. These 
are the ions that were measured by Sembach et al. (1999, 
2000), Collins et al. (2004), and Fox et al. (2005) in the 
high- velocity metal ion absorbers (see e.g. Table [1]). 

4. METAL- ABSORPTION COLUMN DENSITIES 

In this section we present computations of the inte- 
grated metal-ion column densities that are produced in 
the conductively evaporating layers. 

The characteristic column density in an interface, 
A/him = "-him-R, is given by (see equation [5]) , 



Ahim = 2 x 



2khim7him 

250 s ^iHIMC 3 1IM 



T 2 

1 _ 1n i7 J 6,HIM _2 

= 1.7 x 10 — - — cm , 



"000.3 



(19) 

where ^him is the mean mass per particle, T^him is the 
HIM temperature in units of 10 6 K, and </>o.3 = 4> s /0.3. 
This equation shows that in general, models with lower 
saturation parameters produce larger columns. The col- 
umn density of a given ionic species is related to this 
characteristic column by an "ionization correction" , that 
depends on the interface parameters. 

For gas in collisional ionization equilibrium, the local 
ion fractions are a function of the temperature only. The 
ionization correction and implied metal-ion column den- 
sities therefore depend only on a (through its control of 
the temperature profile, see Figuref!]) and on the ambient 
temperature, Thim- However, the functional dependence 
on cto an d ?him may be complicated. 

For non-equilibrium conditions, departures from equi- 
librium ionization and the ion fractions at a given tem- 
perature depend on the ratio of the heating time-scale 
and the ionization time-scale. Equation [17] shows that 
in the absence of photoionization, the non-equilibrium 
abundances and columns still depend only on cto and 
Thim- The inclusion of photoionization by the meta- 
galactic radiation field introduces an implicit depen- 
dence on the gas density (or ionization parameter, see 
equation I18[) . The non-equilibrium column densities are 
therefore general functions of all three interface paramo- 
terQ 

Here we wish to find conditions that yield metal-ion 
column densities comparable to those observed in the 

14 We find that the column densities listed in Table [4] cannot be 
expressed as simple power-law fits in Tjjimi P/^b, and -R. 
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Fig. 7. — (a) Non-equilibrium ionization in a high-pressure model with R = 1.3 kpc, Phtm/^ = 50 cm -3 K, and Thim = 2 x 10 6 K. The 
solid and dashed lines show the non-equilibrium and equilibrium abundances, respectively. The non-equilibrium abundances reach their 
peaks further away from the cloud center and span larger path-lengths, (b) The effect of photoionization by an external radiation field, for 
low-pressure model with R = 7.5 kpc, Phim/^ = 0.1 cm -3 K, and Phim = 2 x 10 K. The solid and dashed lines show the abundances 
including photoionization by the metagalactic radiation file, and in the absence of external radiation, respectively. Photoionization increases 
the level of ionization in the interface. 



TABLE 3 
Model Parameters 



Phim 


R 


P/k B 


U* 




m cl 




max(M)t 




(K) 


(kpc) 


(cm- 3 K) 






(gs- 1 ) 




Gyr 


10 6 


3 


0.1 


9.7 x 10~ 2 


85.62 


6.45 x 10 22 


0.03 


0.60 


330 


10 6 


3 


1 


9.7 x 10~ 3 


8.83 


6.64 x 10 22 


0.20 


0.58 


47.8 


10 6 


3 


5 


1.9 x 10~ 3 


1.81 


6.80 x 10 22 


0.63 


0.38 


14.8 


10 6 


3 


10 


9.7 x 10~ 4 


0.91 


6.86 x 10 22 


0.80 


0.23 


11.6 


10 6 


3 


20 


4.8 x 10~ 4 


0.46 


6.93 x 10 22 


0.89 


0.13 


10.3 


10 6 


3 


50 


1.9 x 10~ 4 


0.19 


7.02 x 10 22 


0.96 


0.06 


9.4 


2 x 10 6 


1.3 


50 


1.9 x 10~ 4 


3.34 


1.67 x 10 23 


0.43 


0.50 


8.8 


2 x 10 6 


3 


0.1 


9.7 x 10~ 2 


666.8 


3.56 x 10 23 


0.01 


0.60 


178 


2 x 10 6 


3 


1 


9.7 x 10~ 3 


68.70 


3.66 x 10 23 


0.04 


0.60 


43.4 


2 x 10 6 


3 


5 


1.9 x 10~ 3 


14.04 


3.74 x 10 23 


0.13 


0.59 


13.1 


2 x 10 6 


3 


10 


9.7 x 10~ 4 


7.08 


3.77 x 10 23 


0.24 


0.57 


7.0 


2 x 10 6 


3 


20 


4.8 x 10~ 4 


3.58 


3.81 x 10 23 


0.41 


0.51 


4.1 


2 x 10 6 


3 


50 


1.9 x 10~ 4 


1.45 


3.86 x 10 23 


0.69 


0.33 


2.4 


2 x 10 6 


4.2 


1 


9.7 x 10~ 3 


49.40 


5.09 x 10 23 


0.05 


0.60 


24.9 


2 x 10 6 


7.5 


0.1 


9.7 x 10~ 2 


266.7 


8.87 x 10 23 


0.01 


0.60 


71.6 


2 x 10 6 


7.5 


1 


9.7 x 10~ 3 


27.48 


9.14 x 10 23 


0.08 


0.60 


8.7 


2 x 10 6 


7.5 


5 


1.9 x 10~ 3 


5.62 


9.34 x 10 23 


0.29 


0.56 


2.3 


2 x 10 6 


7.5 


10 


9.7 x 10~ 4 


2.83 


9.43 x 10 23 


0.49 


0.47 


1.4 


2 x 10 6 


7.5 


20 


4.8 x 10~ 4 


1.43 


9.52 x 10 23 


0.70 


0.33 


0.95 


2 x 10 6 


7.5 


50 


1.9 x 10~ 4 


0.58 


9.64 x 10 23 


0.87 


0.16 


0.77 


10 7 


3 


10 


9.7 x 10~ 4 


831.1 


1.98 x 10 25 


0.007 


0.60 


4.7 


10 7 


3 


20 


4.8 x 10~ 4 


419.3 


2.00 x 10 25 


0.008 


0.60 


4.0 


10 7 


3 


50 


1.9 x 10~ 4 


169.7 


2.02 x 10 25 


0.02 


0.60 


1.9 


10 7 


7.5 


10 


9.7 x 10~ 4 


332.4 


4.95 x 10 25 


0.01 


0.60 


1.3 


10 7 


7.5 


20 


4.8 x 10~ 4 


167.7 


4.99 x 10 25 


0.02 


0.60 


0.75 


10 7 


7.5 


50 


1.9 x 10~ 4 


67.9 


5.05 x 10 25 


0.04 


0.60 


0.35 



Note. — 
Results are shown for <j>s = 0.3. 

* - The ionization parameter U = 9.7 x 10~ 3 (P/Aib) -1 - 
f - The maximal Mach number in the interface, max(M). 

| - The evaporation time, t eV ap = mr/rh Gyr, where m-j is the warm cloud gas mass in units of 10 7 Mq, and rh = umi c \. 
The highlighted lines correspond to the models considered in sections 4.1 and 4.2 below. 
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Fig. 8.1. — Ion fractions as a function of scaled radius r/R — 1, 
in conductive interfaces surrounding dwarf-galaxy-halos. The am- 
bient temperature, Thim; cloud radius, R, and bounding pressure, 
P are indicated in each panel. In Fig 8.1 the non-equilibrium 
ion fractions are for Tjjim = 2 x 10 K, R = 1.3 kpc, and 
P/k B = 50 cm~ 3 K. 
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Fig. 8.2.— Same as figure 8.1 but for T mM = 2 X 10 6 K, R = 
4.2 kpc, and P/k B = 1 cm" 3 K. 

high-velocity ionized absorbers, including the high ob- 
served O VI column (iVovi ~ 10 14 cm -2 ). We focus on 
the column densities of the high ions C IV, N V, O VI, 
Si III, Si IV, and S III, as a function of the impact pa- 
rameter from the cloud center, b. In Table [4] we list the 
(two-sided) central (b = 0) and maximal column densi- 
ties for the different models considered in section 3. For 
the evaporating gas, the line-of-sight column density of 
ion i of element m observed at an impact parameter b is, 

Ni = 2 r n u (r)A mXl (r)rdr ^ 

imax(fc,i?) Vr 2 -b 2 

where A m is the abundance of element m relative to hy- 
drogen, nn(r) is the position dependent hydrogen den- 
sity, and Xi(r) are the position-dependent ion fractions, 
as presented in section[3] The maximal columns are eval- 
uated for each ion independently, and occur at different 
impact parameters for each ion and model. 

The results presented in Table 0] are computed for Z = 
0.1 times solar. Since the fractional ion abundances are 
independent of metallicity, the integrated columns are 
proportional to Z. (However, as discussed in 12.1.21 the 



range of Co values for which a self-consistent evaporating 
solution can be found is a function of Z .) 

Table [4] lists the central column densities produced in 
the interface layers. When observing evaporating clouds, 
the observed absorption line properties will generally de- 
pend not only on the conductive interface, but also on 
the ionization properties of the embedded photoionized 
clouds. In the text that follows we discuss, as illustrative 
examples, the total (cloud + interface) column densi- 
ties that arise in two specific models. First we consider a 
low-ionization cloud embedded in a high-pressure corona, 
with P/k-Q = 50 cm -3 K (§4.1), and then we study 
a more highly-ionized cloud embedded a low-pressure 
medium, with P/kB = 1 cm -3 K (in §4.2). The lines 
representing these two illustrative models are highlighted 
in Table El 

4.1. A Low- Ionization Cloud Embedded in a High 
Pressure Medium 

As an example for a low-ionization photoionized cloud 
embedded in a high-pressure corona, we rely on the phys- 
ical properties of the "CHVC-scale" model presented in 
GS04 ("model A" in GS04). This model represents a 
10 s M Q dark-matter halo that contains a warm gas mass 
of 1.1 x 10 6 M , which is subjected to an external bound- 
ing pressure, P/k-Q = 50 cm -3 K. The resulting pho- 
toionized cloud has a radius R = 1.3 kpc. The ionization 
parameter on the cloud surface is low, U = 1.9 x 10~ 4 , 
and high ions are thus inefficiently produces in this cloud. 
Here we assume that the cloud is surrounded by a hot 
Galactic corona, with a temperature Thim = 2 x 10 6 K, 
and photoionized by the metagalactic radiation field. 

These parameters imply that the cloud evaporates in a 
marginally-saturated flow, with Co ~ 3.3. The mass loss 
rate is suppressed by a factor 0.43 relative to that in a 
classical diffusive flow, and the maximal Mach number 
in the flow is 0.50. 

Figure 8.1 displays the ion fractions that arise in the 
conductive interface in this model. In evaluating the 
absorption-line signatures that arise from the evaporat- 
ing cloud, we use the abundances given in GS04 for 
"model A" (Figure 3 in GS04) for the ion fractions within 
the photoionized cloud, and the ion fractions computed 
here and presented in Figure 8.1 for the evaporating lay- 
ers. We then evaluate the projected column densities as 
a function of impact parameter from the cloud center. 
These column densities include the contributions of both 
the photoionized cloud and the conductive interface sur- 
rounding it. 

Figure [9] shows the projected column densities in our 
evaporating high-pressure model. The ionization in this 
conductive interface results in a maximal O VI column 
density of 3.2 x 10 12 cm -2 (compared to 4.5 x 10 6 cm -2 
in the purely photoionized cloud), and a C IV column 
density of 2.5 x 10 12 cm -2 (compared to 2.6 x 10 11 cm -2 
in the purely photoionized cloud). The evaporating lay- 
ers increase the O VI column by ~ 6 orders of mag- 
nitude, and the C IV column by ~ 1 order of magni- 
tude relative to the purely photoionized cloud discussed 
in GS04. However, for this model, the enhanced columns 
still do not bring the model into harmony with the ob- 
servations of the high-velocity metal-ion absorbers (see 
e.g. Table [J). 
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Fig. 9. — Projected column densities as a function of impact 
parameter in a high-pressure model. For this model Ft = 1.3 kpc, 
P/k B = 50 cm" 3 K, and T mM = 2 X 10 6 K. 



For pressure supported photoionized clouds in high- 
pressure environments (P ~ 50 cm~ 3 A"), the ionization 
parameter is generally too low to allow for efficient pro- 
duction of high ions such as C 3+ and 5+ . The high-ion 
column densities are therefore dominated by ions pro- 
duced in the evaporating layers, and are not significantly 
affected by the details of the cloud's ionization state. 

4.2. A Highly-ionized Cloud Embedded in a Low 
Pressure Medium 

As an illustrative example for the physical proper- 
ties of a more highly ionized photoionized cloud em- 
bedded in a low-pressure medium, we rely on the stan- 
dard dwarf-scale model presented in GS04 ("model C" 
in GS04), based on the observed properties of Local 
Group dwarf galaxies Leo A and Sag DIG (Sternberg et 
al. 2002). This model represents a 2 x 10 9 M Q dark- 
matter halo that contains a total warm gas mass of 
1.6 x 10 7 M , and is subjected to an external bound- 
ing pressure P/fce = 1 cm~ 3 K. The resulting cloud has 
a radius R <~ 4.2 kpc. The ionization parameter on the 
surface of this cloud is higher, U = 9.7x 10 -3 , and high 
ions such as C 3+ and N 4+ are efficiently produced in the 
ionized envelopes of this cloud. Here we assume that the 
halo is surrounded by a hot intergalactic medium, with 
a temperature Thim = 2 x 10 6 K, and is photoionized by 
the metagalactic radiation field. 

This cloud then evaporates in a saturated flow, with 
Co ~ 50. The mass- loss rate is suppressed by a factor 
0.05 relative to a classical diffusive flow, and the maximal 
Mach number in the flow is 0.60. 

Figure 8.2 displays the ion fractions of the low ions 
C II, N I, O I, Si II, and S II, and the high-ions C IV, 
N V, O VI, Si III, Si IV, and S III in the conductive 
interface surrounding this dwarf scale halo. 

For the column densities produced in this model, we 
rely on the ion-abundances given in GS04 ("model C", 
Figure 5) for the photoionized component, and the ion 
fractions computed here and presented in Figure 8.2 for 
the conductive interface surrounding it. We compute the 
projected column densities as a function of impact pa- 
rameter from the cloud center, again including both the 
photoionized core and the surrounding evaporating en- 
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Fig. 10. — Projected column densities as a function of impact 
parameter in a dwarf galaxy-scale model. For this model R = 
4.2 kpc, P/k B = 1 cm" 3 K, and T mM = 2 X 10 6 K. 

Figure [TU] shows the projected column densities in our 
dwarf-galaxy-scale model. The ionization in this conduc- 
tive interface results in a maximal O VI column density 
of 5.5 x 10 11 cm~ 2 (compared to 2.9 x 10 11 cm" 2 in the 
photoionized cloud), and a maximal C IV column den- 
sity of 8.3 x 10 12 cm~ 2 , comparable to that produced in 
the photoionized cloud). This is only a minor enhance- 
ment relative to the purely photoionized cloud, and the 
columns in this dwarf galaxy-scale model are therefore 
still too low to account for the observed properties of the 
metal- ion absorbers (e.g. Table [TJ. 

4.3. Diagnostics 

Diagnostic diagrams for evaporating gas may be con- 
structed using Equation [20] and the computational data 
presented in § [3j and used to probe the ionization mech- 
anisms at play in observed metal-ion absorbers. As 
an example, in Figure [TT] we display trajectories for 
Ac iv /No vi versus An v/Ao vi ; as a function of im- 
pact parameter, b through the evaporative interface. The 
impact parameter is represented by the color along the 
curves, from the cloud center where 6 = (blue) to the 
surrounding HIM where b > R (red). The different pan- 
els are for non-equilibrium evaporation for a range of 
dwarf-scale objects, drawn from the range of cloud radii, 
bounding pressures, and HIM temperatures considered in 
Table 01 These diagrams show the possible range of ion 
ratios produced in evaporating dwarf-scale models. Since 
the column densities are proportional to the metallicity, 
the ion ratios are independent of Z. Note that these 
ion ratios do not include any contribution from the em- 
bedded warm cloud, which depends on the astrophysical 
setting. 

5+ is the most highly ionized of the three ions that 
we consider. Nq iv/Ao vi an d N^ v/No vi are there- 
fore large at low temperatures where the 5+ fraction 
is small, and decrease at higher temperature as 5+ 
becomes more abundant. Small impact parameters (in 
blue) sample the lower-T parts of the interface, where 
Ac rv /No vi and An v/Ao vi are large. The cool- 
ing trajectories therefore start in the upper right cor- 
ner of the parameter space. Larger impact parame- 
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Note. — Assuming <f> s = 0.3. Columns were computed for Z = 0.1 times solar, and are proportional to Z. 



ters only sample the hotter parts of the interface, and 
the trajectories evolve in the direction of decreasing 
Ac rv/Ao vi and An v/Ao vi- Finally, as the temper- 
ature approaches Trim (2 x 10 6 K in Figure [H]), 5+ is 
replaced by higher ionization states, and the evolutionary 
tracks turn up again. 

An additional ion that is commonly observed in high- 
velocity absorbers, and that has been used as a diagnostic 
for the physical conditions prevailing in these systems, is 
Si III (Collins et al. 2004; Collins et al. 2009; Shull et 



al. 2009). The typical observed column densities of Si III 
in the high-velocity absorbers are ~ 3 x 10 13 cm~ 2 (Shull 
et al. 2009). Si III is never efficiently produced in evapo- 
rating dwarf-galaxy-scale clouds. For low bounding pres- 
sures the ionization parameter is too high to efficiently 
produce Si III in the evaporating layers, but even at high 
pressures (P ~ 50 cm -3 K), for which Si III is the dom- 
inant ionization state at the cloud surface, columns in 
excess of 10 12 cm -2 are rarely produced. Si III must 
therefore be produced by alternative ionization mecha- 
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Fig. 11. — Column density ratios Nq iv/Nq VI versus v/No vi f° r evaporating gas. The impact parameter b is indicated by the 
color along the trajectories, from the cloud center (blue) to the surrounding HIM (red). The different panels are for different interface 
parameters, as indicated by the legends. These ion ratios do not include any contribution from the embedded warm cloud. 



nisms. It may be efficiently produced by photoionization 
within an embedded warm cloud, as we demonstrate in 
Figures E] and HO] 

4.4. High O VI columns in evaporating clouds? 




10 



10 



10 



Fig. 12. — Central (two-sided) O VI column density versus sat- 
uration. Models inside the shaded areas are not self-consistent, as 
cooling may not be neglected (see section ^. 1.211 . 

The results presented in Table Q] and illustrated by the 
examples discussed above, show that the typical O VI 
columns that are produced in the different dwarf-scale 
models are < 10 13 cm -2 . Can column densities of or- 
der 10 14 cm -2 be produced in conductively evaporating 
clouds? 



In Figure [T2l we show the central (two-sided) O VI col- 
umn density as a function of saturation parameter for 
our full grid of models, assuming Z = 0.1. The range of 
saturation parameters for which no self-consistent evap- 
orating solution can be found has been shaded in the 
plot. In the shaded region, cooling overcomes conduc- 
tive heating, evaporation stops, and instead the ambient 
medium condenses onto the clouds. Our theory does not 
apply to condensing clouds. Figure [12] shows that cen- 
tral O VI column densities of order 10 14 cm -2 cannot 
be obtained for any Z = 0.1 evaporating model. The 
largest central O VI columns, obtained for cro ~ 0.2, are 
- 2 x 10 12 cm" 2 . 

We note that for solar metallicity clouds, the O VI 
columns are a factor of 10 larger. However, since for Z = 
1 self-consistent evaporating solutions can only be found 
for er > 0.5, central O VI columns of order 10 14 cm~ 2 
can still not be obtained. 

Since the observed column density is a function of the 
impact parameter, we also computed the maximal O VI 
column density obtained in each self-consistent evapo- 
rating model. Inspecting our entire parameter space 
(0.1 < P/k B < 10 4 cm" 3 K, 1 pc < R < 100 kpc, 
5 x 10 5 < T < 10 7 K), we find that for Z = 0.1, max- 
imal column densities of order 10 14 cm -2 only occur in 
evaporating clouds (cro > 0.2) for extreme combinations 
of parameters, with T ~ 10 7 K, and P ~ 10 4 cm -3 K. 
We therefore conclude that non-equilibrium ionization in 
conductively evaporating clouds an unlikely explanation 
of the high-ions observed in the high-velocity metal-ion 
absorbers. An alternative explanation is that O VI is 
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created in a flow of cooling gas (e.g. Edgar & Cheva- 
lier 1986; Heckman et al. 2002; Gnat & Sternberg 2007; 
Gnat & Sternberg 2009). 

While a conductive origin for the O VI absorption ob- 
served in the ionized high-velocity clouds appears un- 
likely, recent observations indicate that evaporating O VI 
may be present in more local environments. Savage & 
Lehner (2006) report on FUSE absorption line obser- 
vations toward a sample of local white dwarfs. O VI 
was detected along 24 lines of sight, with column den- 
sities between 2.5 x 10 12 and 4.0 x 10 13 cm" 2 . Based 
on the observed columns and on the kinematics of the 
absorption, Savage & Lehner suggest that this sample of 
local absorbers may probe conductive interfaces in the 
local interstellar medium. Their analysis relied on the 
plane-parallel conduction-front models of Borkowski et 
al. (1990). Our spherical models for typical interstel- 
lar parameters (T ~ 10 6 K, P ~ 4 x 10 3 cm" 3 K, and 
R ~ 15 pc) support this interpretation. For interstellar 
conditions the saturation parameter <7o ~ 0.5, and we 
find that O VI column densities of order 10 13 cm" 2 are 
expected for solar metallicity gas. The predicted Doppler 
widths (15—30 km s" 1 ) and flow velocities (5—55 km s" 1 ) 
are also consistent with the Savage & Lehner observa- 
tions. 

5. SUMMARY 

In this paper we present computations of the non- 
equilibrium ionization states of H, He, C, N, O, Si, and 
S in thermally conductive interfaces, surrounding warm 
(and photoionized) gas clouds that undergo steady evap- 
oration into a hot ambient medium. 

We rely on the formalism developed by Dalton & Bal- 
bus (1993; DB93) for solving the temperature profile in 
the evaporating gas. This formalism allows the thermal 
conductivity to change continuously from the classical 
diffusive form (Spitzer 1962) to a flux-limited saturated 
form (Cowie & McKee 1977). The temperature profile 
then depends only on the global level of saturation in the 
flow, as described by the saturation parameter, (Jo- 
in their analytic solution for the temperature profile, 
DB93 neglect the kinetic energy term as well as the ra- 
diative losses term in the energy equation. As opposed 
to previous claims, we find that the neglect of the kinetic 
energy term is accurate only for classical, unsaturated 
flows, and applies to saturated flows only if the parame- 
ter <f) a (see Equation 2) is less than 0.6. 

We explicitly solve the time-dependent ionization 
equations for the evaporating gas, including also pho- 
toionization by the present-day metagalactic field. As 
the gas evaporates, heating can become fast compared 
to the ionization time-scale, and the gas tends to remain 
under-ionized. We study how departures from equilib- 
rium ionization affects the ion distributions and result- 
ing cooling rates in the evaporating layers. We show 
that the existence of underionized species, most signifi- 
cantly He + , significantly enhances the cooling rates. We 
estimate numerically the range of saturation parameters 
for which radiative cooling plays a significant role in de- 
termining the dynamics of the gas. We find that for a 
0.1 (1) times solar metallicity gas, radiative losses are 
significant for <7q i$ 0.2 (0.5). For larger values of cto 
the importance of radiative cooling is generally smaller, 
and depends on the neutral fraction at the surface of the 



evaporating cloud. 

We use our conduction models to extend the purely 
photoionized cloud models for local high-velocity metal 
line absorbers that we presented in Gnat & Stern- 
berg (2004; GS04). In that paper, we modeled the ab- 
sorbers as warm (~ 10 4 K) clouds that are photoion- 
ized by the metagalactic radiation field, and embedded 
in low-mass (< 2xl0 9 M Q ) dark-matter mini-haos. Here, 
we consider the additional ionization that occurs in the 
conductive envelopes, as these clouds evaporate into a 
hot corona and/or intergalactic medium. We consider a 
wide range of input parameters (cloud radius, thermal 
pressure and temperature of the hot medium) as sum- 
marized in Table [3l This table also lists the ionization 
parameter, global saturation parameter, mass-loss rate, 
maximal Mach number, and evaporation time scale for 
each model. We assume that <fi s = 0.3, so that the flow 
is everywhere subsonic. 

In §3 we compute the ionization states as a function 
of distance from the evaporating cloud surface. We plot 
the resulting ion fractions in Figure 8. In §4 we present 
the metal absorption column densities arising in conduc- 
tively evaporating minihalo clouds. We list the column 
densities of the high ions that are produced in the in- 
terface layers of our dwarf-scale models in Table |4] We 
then focus on two illustrative models, for which we in- 
clude the contributions to the ion distributions from both 
a central photoionized cloud and the conductive inter- 
face surrounding it. In §4.1 we discuss the ion fractions 
and projected column densities in a low-ionization model 
embedded in a high-pressure environment. In §4.2, we 
present results for a more highly ionized cloud embedded 
in a low-pressure medium. We find that the contribu- 
tion of the conductive interface enhances the formation 
of high-ions such as C IV and O VI. However, in all cases 
the enhanced columns are still too low to account for 
the O VI columns observed in the high-velocity metal- 
ion absorbers (e.g. Tabled]). The central O VI column 
density for evaporating clouds is limited to < 10 13 cm" 2 . 
These limits are reached for small saturation parameters, 
just at the point where radiative losses become significant 
enough to result in condensation. While a conductive ori- 
gin for the O VI absorption observed in the ionized high- 
velocity clouds appears unlikely, our models do support 
the conclusion by Savage & Lehner (2006) that evapo- 
rating O VI absorption occurs in the local ISM, with 
characteristic columns of ~ 10 13 cm -2 . 

We use the predictions for the high-ion column densi- 
ties produced in evaporating clouds (see Table [3} to con- 
struct ion-ratio diagnostic diagrams. In Section 4.3, we 
discuss one example, Nq iv/Ao vi versus An v/Ao vij 
and show how these ratios vary with impact parameter 
in the evaporating layers, and how they may be used 
to probe the ionization processes affecting metal-line ab- 
sorbers. 
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APPENDIX 

ON THE SENSITIVITY TO THE SPECTRAL SLOPE NEAR THE LYMAN LIMIT 

Our standard fit for the metagalactic radiation field is given by Equation [13J and displayed in Figure 0] (solid dark 
curve). Near the Lyman limit, this representation varies as i/~ 313 . However, the spectral slope may be significantly 
shallower at the Lyman limit, with J„ varying as u , and steepen gradually with increasing frequency (e.g. Shull 
et al. 1999; Telfer et al. 2002; Scott et al. 2004; Shull et al. 2004; Zheng et al. 2004; Faucher-Giguere et al. 2009). In 
this Appendix we present results for the ionization properties in evaporating conduction fronts for a "maximal field" 
(Gnat & Sternberg 2004) that is represented by, 



Jv 
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2.512 x W- 2 J u0 (% 
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< 18.4 



(Al) 



Hz 1 sr , and 



and displayed by the dashed curve in FigureHJ As for our standard field, J„q = 2 x 10~ 23 erg s _1 cr 
vq is the Lyman limit frequency. Our standard and maximal fields span the range of possible spectral slopes suggested 
by observations and theory. 

We repeat the calculations presented in Sections (HO andHJassuming irradiation by the maximal field. The increased 
intensity above the Lyman limit yields higher photoionization rates, which affect the column densities of intermediate 
ions. In Table [5j we list the (two-sided) central (b = 0) and maximum column densities for the different models 
considered in §3, assuming photoionization by the maximal field. A comparison of Tables U and [S] demonstrates that 
the high-ion column densities are generally insensitive to the exact value of the spectral slope near the Lyman limit. 

Table [S] confirms that low pressure/density models are more affected by the increased UV intensity than higher 
pressure models, and that intermediate ions are more affected than the high ions. In all cases, we find that the 
columns computed with the standard field agree with those computed with the maximal field to better than 0.5 dex. 
As we describe below, in most cases the agreement is better. 




P (cm" J K) 

Fig. 13. — The range of A r max /A'standard ratios as a function of the external bounding pressure, for Si III, C IV, N V, and O VI produced 
in dwarf-scalc-halos with T = 1 — 2 X 10 6 K. Here -/V stan< j ar d are the columns computed with the standard field described in Section l2.2l and 
JVmax are the columns computed with the maximal field, which is shallower above the Lyman limit. The high-ion columns are insensitive 
to spectral slope near the Lyman limit, and the intermediate-ions are only moderately affected, mostly for low bounding pressures. 

Figure [13] compares the columns computed with the maximal field, N ma-X , to the column computed with the standard 
field, iVstandard- We display the range of -/V max /Ar st andard ratios as a function of the external bounding pressure, for 
Si III, C IV, N V, and O VI produced in dwarf scale halos with T = 1 - 2 x 10 6 K. 

Figure [13] shows that for P > 1 cm~ 3 K, the O VI and N V columns remain similar, the C IV columns agree to 
within 20%, and the Si III columns agree to within 30%. For the lowest pressure that we consider, P = 0.1 cm~ 3 K, 
the differences are larger, but are still within a factor 2.5. We conclude that column densities of the high ions are only 
weakly dependent on the exact spectral slope near the Lyman limit. 
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TABLE 5 

Column Densities for the Maximal Radiation Field 





R 


P/k B 






C IV 


N V 


O VI 


Si III 


Si IV 


S III 


10 6 


3 


0.1 


85.62 


central 


6.2 x 


10 9 


6.2 


X 


10 9 


8.3 x 


10 10 


1.5 


X 


10 5 


6.4 


X 


10 6 


8.0 x 10 4 










maximal 


2.6 x 


10 10 


1.8 


x 


10 1() 


1.9 x 


10 11 


1.0 


X 


10 6 


4.1 


X 


10 7 


3.6 x 10 s 


10 6 


3 


1 


8.83 


central 


1.4 x 


10 11 


4.9 


x 


10 10 


3.9 x 


10 11 


4.1 


X 


10 8 


3.3 


X 


10 9 


2.0 x 10 8 










maximal 


5.3 x 


10 11 


8.5 


x 


10 1() 


5.6 x 


10 11 


5.4 


x 


10 9 


2.7 


x 


10 10 


2.3 x 10 9 


10 6 


3 


5 


1.81 


central 


3.4 x 


10 11 


1.3 


x 


10 11 


1.1 X 


10 12 


5.0 


x 


10 9 


1.9 


x 


10 1() 


2.0 x 10 9 










maximal 


1.2 x 


10 12 


2.8 


X 


10 11 


2.0 x 


10 12 


9.6 


X 


10 10 


1.5 


X 


10 11 


3.4 x 10 10 


10 6 


3 


10 


0.91 


central 


4.3 x 


10 11 


1.7 


X 


10 11 


1.7 x 


10 12 


1.1 


X 


10 10 


2.7 


X 


10 1() 


5.1 x 10 9 










maximal 


2.1 x 


10 12 


5.0 


x 


10 11 


3.6 x 


10 12 


2.8 


x 


10 11 


2.5 


x 


10 11 


1.2 x 10 11 


10 6 


3 


20 


0.46 


central 


5.2 x 


10 11 


2.2 


X 


10 11 


2.5 x 


10 12 


1.2 


X 


10 10 


3.3 


X 


10 10 


6.5 x 10 9 










maximal 


3.6 x 


10 12 


8.5 


x 


10 11 


6.2 x 


10 12 


3.6 


X 


10 11 


4.4 


X 


10 11 


1.7 x 10 11 


10 6 


3 


50 


0.19 


central 


7.1 x 


10 11 


3.4 


x 


10 11 


4.5 x 


10 12 


5.3 


X 


10 10 


4.6 


X 


10 1() 


3.0 x 10 10 










maximal 


6.7 x 


10 12 


1.6 


X 


10 12 


1.2 x 


10 13 


1.8 


X 


10 12 


7.9 


X 


10 11 


1.1 x 10 12 


2 x 10 6 


1.3 


50 


3.34 


central 


4.1 x 


10 11 


1.6 


X 


10 11 


1.3 x 


10 12 


7.5 


X 


10 9 


2.5 


X 


10 1() 


3.6 x 10 9 










maximal 


2.5 x 


10 12 


5.4 


X 


10 11 


3.2 x 


10 12 


2.3 


X 


10 11 


3.3 


X 


10 11 


1.0 x 10 11 


2 x 10 6 


3 


0.1 


666.8 


central 


5.4 x 


10 9 


5.6 


X 


10 9 


6.9 x 


10 10 


1.3 


X 


10 5 


4.9 


X 


10 6 


7.1 x 10 4 










maximal 


2.6 x 


10 10 


1.9 


x 


10 1() 


1.9 x 


10 11 


9.8 


X 


10 5 


3.5 


X 


10 7 


5.3 x 10 s 


2 x 10 6 


3 


1 


68.70 


central 


1.3 x 


10 11 


3.7 


X 


10 1() 


2.4 x 


10 11 


3.8 


X 


10 8 


3.3 


X 


10 9 


1.6 x 10 s 










maximal 


4.4 x 


10 11 


6.5 


X 


10 10 


3.2 x 


10 11 


4.3 


X 


10 9 


2.4 


X 


10 1() 


1.5 x 10 9 


2 x 10 6 


3 


5 


14.04 


central 


3.0 x 


10 11 


9.0 


X 


10 10 


6.2 x 


10 11 


4.6 


X 


10 9 


1.7 


X 


10 1() 


1.7 x 10 9 










maximal 


9.6 x 


10 11 


1.8 


X 


10 11 


9.5 x 


10 11 


9.0 


X 


10 10 


1.4 


X 


10 11 


3.1 x 10 1() 


2 x 10 s 


3 


10 


7.08 


central 


3.6 x 


10 11 


1.2 


X 


10 11 


9.0 x 


10 11 


6.1 


X 


10 9 


2.1 


X 


10 1() 


2.7 x 10 9 










maximal 


1.5 x 


10 12 


3.0 


x 


10 11 


1.6 x 


10 12 


1.5 


X 


10 11 


1.9 


X 


10 11 


6.4 x 10 1() 


2 x 10 6 


3 


20 


3.58 


central 


4.1 x 


10 11 


1.5 


X 


10 11 


1.3 x 


10 12 


1.6 


X 


10 10 


2.6 


X 


10 1() 


7.7 x 10 9 










maximal 


2.4 x 


10 12 


5.1 


X 


10 11 


2.9 x 


10 12 


6.6 


X 


10 11 


3.2 


X 


10 11 


3.2 x 10 11 


2 x 10 6 


3 


50 


1.45 


central 


4.7 x 


10 11 


1.9 


x 


10 11 


2.0 x 


10 12 


8.2 


X 


10 9 


2.9 


X 


10 1() 


3.9 x 10 9 










maximal 


4.7 x 


10 12 


1.1 


X 


10 12 


6.8 x 


10 12 


3.6 


X 


10 11 


6.0 


X 


10 11 


1.6 x 10 11 


2 x 10 6 


7.5 


0.1 


226.7 


central 


7.7 x 


10 9 


8.6 


x 


10 9 


1.1 X 


10 11 


1.8 


X 


10 5 


6.7 


X 


10 6 


1.0 x 10 s 










maximal 


4.5 x 


10 10 


3.2 


X 


10 1() 


3.2 x 


10 11 


1.8 


X 


10 6 


6.7 


X 


10 7 


6.6 x 10 5 


2 x 10 6 


7.5 


1 


27.48 


central 


1.8 x 


10 11 


6.3 


X 


10 10 


4.4 x 


10 11 


4.5 


X 


10 8 


3.4 


X 


10 9 


2.4 x 10 s 










maximal 


9.5 x 


10 11 


1.5 


X 


10 11 


6.2 x 


10 11 


8.8 


X 


10 9 


4.2 


X 


10 10 


4.1 x 10 9 


2 x 10 6 


7.5 


5 


5.62 


central 


3.7 x 


10 11 


1.3 


X 


10 11 


1.0 x 


10 12 


6.3 


X 


10 9 


1.9 


X 


10 1() 


2.8 x 10 9 










maximal 


2.4 x 


10 12 


3.7 


X 


10 11 


2.0 x 


10 12 


2.4 


X 


10 11 


2.9 


X 


10 11 


1.1 x 10 11 


2 x 10 s 


7.5 


10 


2.83 


central 


4.0 x 


10 11 


1.6 


x 


10 11 


1.4 x 


10 12 


4.1 


X 


10 9 


2.1 


X 


10 1() 


1.9 x 10 9 










maximal 


2.8 x 


10 12 


6.0 


X 


10 11 


3.6 x 


10 12 


1.6 


X 


10 11 


3.2 


X 


10 11 


6.9 x 10 1() 


2 x 10 6 


7.5 


20 


1.43 


central 


4.7 x 


10 11 


1.9 


X 


10 11 


1.9 x 


10 12 


1.9 


X 


10 10 


3.1 


X 


10 10 


9.2 x 10 9 










maximal 


4.7 x 


10 12 


1.1 


X 


10 12 


6.8 x 


10 12 


1.1 


X 


10 12 


5.9 


X 


10 11 


5.5 x 10 11 


2 x 10 6 


7.5 


50 


0.58 


central 


5.6 x 


10 11 


2.4 


X 


10 11 


2.9 x 


10 12 


3.6 


X 


10 10 


4.0 


X 


10 1() 


1.8 x 10 10 










maximal 


9.5 x 


10 12 


2.1 


x 


10 12 


1.5 x 


10 13 


2.4 


X 


10 12 


1.2 


X 


10 12 


1.3 x 10 12 


10 7 


3 


10 


831.1 


central 


3.7 x 


10 11 


1.3 


X 


10 11 


9.8 x 


10 11 


7.0 


X 


10 9 


1.7 


X 


10 1() 


2.5 x 10 9 










maximal 


2.1 x 


10 12 


3.5 


X 


10 11 


1.8 x 


10 12 


2.2 


X 


10 11 


1.9 


X 


10 11 


7.5 x 10 10 


10 7 


3 


20 


419.3 


central 


4.0 x 


10 11 


1.5 


X 


10 11 


1.2 x 


10 12 


1.2 


X 


10 10 


2.4 


X 


10 1() 


6.0 x 10 9 










maximal 


2.4 x 


10 12 


4.9 


X 


10 11 


2.6 x 


10 12 


6.1 


X 


10 11 


3.1 


X 


10 11 


3.0 x 10 11 


10 7 


3 


50 


169.7 


central 


4.3 x 


10 11 


1.7 


X 


10 11 


1.5 x 


10 12 


7.5 


X 


10 9 


2.5 


X 


10 1() 


3.4 x 10 9 










maximal 


4.1 x 


10 12 


8.7 


X 


10 11 


5.0 x 


10 12 


3.2 


X 


10 11 


5.1 


X 


10 11 


1.4 x 10 11 


10 7 


7.5 


10 


332.4 


central 


4.2 x 


10 11 


1.5 


X 


10 11 


1.3 x 


10 12 


9.7 


X 


10 9 


2.3 


X 


10 1() 


4.1 x 10 9 










maximal 


3.2 x 


10 12 


5.7 


X 


10 11 


3.1 x 


10 12 


4.2 


X 


10 11 


3.7 


X 


10 11 


1.7 x 10 11 


10 7 


7.5 


20 


167.7 


central 


4.4 x 


10 11 


1.7 


X 


10 11 


1.5 x 


10 12 


1.5 


X 


10 10 


2.7 


X 


10 1() 


6.7 x 10 9 










maximal 


4.1 x 


10 12 


8.8 


X 


10 11 


5.0 x 


10 12 


8.4 


X 


10 11 


5.2 


X 


10 11 


3.7 x 10 11 


10 7 


7.5 


50 


67.88 


central 


4.6 x 


10 11 


1.8 


X 


10 11 


1.7 x 


10 12 


2.3 


X 


10 10 


3.0 


X 


10 1() 


1.0 x 10 1() 










maximal 


7.0 x 


10 12 


1.5 


X 


10 12 


9.1 x 


10 12 


1.7 


X 


10 12 


9.2 


X 


10 11 


7.4 x 10 11 



Note. — Assuming <j> s = 0.3. Columns were computed for Z = 0.1 times solar, and are proportional to Z. 
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Electronic (only) Figure Set 8 




Fig. 8.1, 8.2. — Ion fractions as a function of scaled radius r/R — 1, in conductive interfaces surrounding dwarf galaxy-scale halos. The 
ambient temperature, Thim> cloud radius, R, and bounding pressure, P are indicated in each panel. 




Fig. 8.3, 8.4. — Ion fractions versus scaled radius in dwarf galaxy-scale models - continued. 
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Fig. 8.5, 8.6. — Ion fractions versus scaled radius in dwarf galaxy-scale models - continued. 




Fig. 8.7, 8.8. — Ion fractions versus scaled radius in dwarf galaxy-scale models - continued. 




Fig. 8.9, 8.10. — Ion fractions versus scaled radius in dwarf galaxy-scale models — continued. 
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Fig. 8.11, 8.12. — Ion fractions versus scaled radius in dwarf galaxy-scale models - continued. 




Fig. 8.13, 8.14. — Ion fractions versus scaled radius in dwarf galaxy-scale models - continued. 




Fig. 8.15, 8.16. — Ion fractions versus scaled radius in dwarf galaxy-scale models - continued. 
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Fig. 8.17, 8.18. — Ion fractions versus scaled radius in dwarf galaxy-scale models - continued. 




Fig. 8.19, 8.20. — Ion fractions versus scaled radius in dwarf galaxy-scale models - continued. 




Fig. 8.21, 8.22. — Ion fractions versus scaled radius in dwarf galaxy-scale models - continued. 
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